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Abstract 

This is the second paper in the series to introduce a graphical method to loop quantum gravity. We employ the graphical 
method as a powerful tool to calculate the actions of the Euclidean Hamiltonian constraint operator and the so-called inverse 
volume operator on spin network states with trivalent vertices. Both of the operators involve the co-triad operator which contains 
holonomies by construction. The non-ambiguous, concise and visual characters of our graphical method ensure the rigour for our 
calculations. Our results indicate some corrections to the existing results in literatures. 

PACS numbers: 04.60.Pp, 04.60.Ds 


1 Introduction 

It is well known that quantum dynamics is a central issue in loop quantum gravity (LQG) (see [1,2] for review arcticles, and [3,4] 
for books). There are two main approaches to the quantum dynamics, based on the canonical and covariant quantization programs 
respectively. In canonical quantization, the quantum dynamics is determined by some quantum Hamiltonian constraint operator. 
In the covariant program the quantum dynamics is to define a reasonable transition amplitude. One expects that the quantum 
dynamics from the two different approaches can make the same physical predictions. Such an expectation has been achived at 
least in 3-dimensional LQG to certain sense [5]. Although some progress has been made for 4-dimensional case in checking the 
consistency between the two approaches [6, 7, 8], the issue is not yet understood up to now. To understand the relation between 
the canonical and covariant quantum dynamics, we not only need a suitable definition of the Hamiltonian constraint operator, 
hut also have to calculate its matrix elements on given quantum states. In the light of the seminal work by Thiemann [9, 10], 
some mathematically well-defined Hamiltonian constraint operators have been constructed in LQG. For instance, a Hamiltonian 
constraint operator alternative to that in [9] was proposed in [11] and its matrix elements acting on a gauge-invariant trivalent 
vertex was derived by the graphical Penrose binor calculus. Later on, the matrix elements was re-derived in [6], and then the 
formula in [6] was corrected by sign factors in [7] using graphical method. Matter coupling is also an important issue in LQG. In 
the case of gravity coupled to a scalar field, the whole Hamiltonian constraint operator was constructed [10, 12]. The matter part 
of the whole Hamiltonian constraint operator usually contains the “inverse volume operator”, which is defined by the co-triad 
operators. In the symmetric model of loop quantum cosmology (LQC) [13], the analogue of the inverse volume operator is 
bounded above. This fact is sometimes thought as a reason for the singularity resolution in LQC. In particular, it is shown in [14] 
that in spatially curved anisotropic models inverse volume effects may becomes important to bind expansion and shear scalars. 
However, it is shown in [15] that the inverse volume operator with certain ordering in full LQG is unbounded on the zero volume 
eigenstates (at a gauge-invariant trivalent vertex). This throws doubt on whether one can generalize the conclusions of LQC to 
LQG. To definitely understand the inverse volume operator in LQG and its relation to the analogues in certain symmetric models, 
it is necessary to calculate in details its action on the quantum states in LQG. There is no doubt that a simple and practical 
calculation method is desirable to further understand both the inverse volume and the Hamiltonian constraint operators. 

Note that the matrix elements of the Hamiltonian constraint operator were calculated in [11] and [6] by the graphical Penrose 
binor method and the modihed Brink’s graphical method respectively, while the inverse volume operator was calculated in [15] 
by the algebraic method. Although the graphical Penrose binor calculus looks simpler and more intuitive than the algebraic 
calculation, the rules of transforming graphs were not shown in [11]. So it is not obvious whether the rules correspond uniquely 
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to the algebraic manipulations of the formula. On the other hand, a detailed derivation of the matrix elements of the Hamiltonian 
constraint operator was not presented in [6]. 

The graphical method developed by Brink in [16] was extended and introduced to LQG (see e.g. [17]) for a simple and 
non-ambiguous calculation method. This method consists of two ingredients, graphical representation and graphical calculation. 
The algebraic formula can be represented by the corresponding graphical formula in an unique and unambiguous way. Then 
one can do calculations following the simple rules of transforming graphs, corresponding uniquely to the algebraic manipulation 
of the formula. The graphical method displayed a powerful efficiency in the derivation of the matrix element of the volume 
operator which involves only the flux operator in [17]. In this paper, we will consider the actions of the Euclidean gravitational 
Hamiltonian constraint operator and the inverse volume operator on spin network states. Both operators depend also on the 
honolomies in addition to fluxes. Our aim is in two folds. One is to show that our graphical method is suitable to calculating 
the actions of different kinds of operators on spin network states. The other is to cross-check the results obtained in literatures, 
on which some important applications are based. Note that. In order to obtain the on-shell anomaly-free quantum constraint 
algebra, one has to employ degenerate triangulation at the co-planar vertices of spin networks in the regularization procedure of 
Thiemann’s Hamiltonian. This problem has been overcome by a new proposed Hamiltonian constrain operator in [18]. 

This paper is divided into four sections. In section 2, we will give a brief review of the construction of Thiemann’s Euclidean 
Hamiltonian constraint operator, and then calculate its action on gauge invariant bivalent spin network states by our graphical 
method. In section 3, we compute the action of the inverse volume operator appeared in the Hamiltonian constraint for gravity 
coupled to a scalar field. The results are summarized in section 4. 


2 The Hamiltonian constraint operator 

To quantize a classical function, we need to first regularize it into a formula represented by the fundamental variables, then to 
replace the variables by their quantum operators and thus obtain a regularized quantum operator. Einally, we remove regulator 
by a suitable limit and obtain the corresponding quantum operator. In this section, we first recall the construction of Thiemann’s 
Hamiltonian constraint operator, and then derive the action of the Euclidean Hamiltonian constraint operator on a spin network 
function over trivalent vertices. 


2.1 Quantization of the Hamiltonian constraint 


The classical Hamiltonian constraint of pure gravity in the connection formulation of general relativity is given by 


H(N)^ 


1 

2k 


f 


d^xN- ' ^ 


Vdet(g') 



2(1 +p^)KlKi] H^N) - 2(1 + l}^)T{N ), 


( 2 . 1 ) 


where k - 2>nG, p is the Barbero-Immirzi parameter, is the curvature of S U(2) connection A],, the sensitized triad E'j is the 
conjugate momentum of A]., K'^ebi is the extrinsic curvature of a spatial hypersurface E in a spacetime with cotriad e*,-, and det(^) 
denotes the determinant of the 3-metric q^b- The function H^(N) is called the Euclidean Hamiltonian constraint. In what follows, 
we focus on the regularization of H^{N). Let us triangulate E into tetrahedra A so that the above integral becomes a sum of 
integrals over A, i.e., We denote the triangulation of E by T(e). The small parameter e indicates the “length” of the 

edges of A. Eor each A, we single out one of its vertices and called it the base-point v(A) of A and denote its three edges outgoing 
from v(A) by i/(A), / = 1,2,3. Taking the limit e ^ 0 corresponds to sharking A to v(A). Let a/j{A) ;= i/(A) o au o iy*(A) 
be the loop based at v(A), where ajj is the edge of A from the endpoint of S/(A) to the endpoint of Sy(A). Then the Euclidean 
Hamiltonian constraint can be written in the form [4, 9] 
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2 iV(v(A))e"^tr(/!„„(A)/z.^(A){/iJ(A)> 


Aer(E) 


=:lim//^(^)(iV), 


( 2 . 2 ) 


where N := sgn(det(e^))N, A^ A^t/c with Tk{k- 1,2,3) the generators of su{2), he = he{A) is the holonomy of a connection A 
along an edge e, V denotes the volume function of E, and in the second step we have used the identities 


^ijk 
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sgn(det(4))e''''''ec = sgn(det(4))e“*‘'^{A', V}6ki, tr(nT;) = —du . 


(2.3) 
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In what follows we will drop the bar over N. Replacing V by V, holonomies by holonomy operators (since the holonomy operator 
acts as a multiplication operator, we also omit the hat for simplification of notation), and the Poisson bracket by 1 / {ih) times the 
commutator, the Euclidean Hamiltonian constraint operator reads 


3itiK^/3 Si., iV(v(A))e^''^tr(/z„„(A)/z,^(A)y/zJ(A)) 


3ifiK^/3 


iV(v(A))i/f, 


(2.4) 


Aer(E) 


The volume operator V actes on a cylindrical function fy over y by [19, 20] 
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(2.5) 


where = V^, qia, ej, ck) = sgn(det(e/( 0 ), ey( 0 ), exiO))), and 

qiJK -^iSijkJeijijJe^ , (2.6) 

here 7],^ is the self-adjoint operator of the right-invariant vector field on the copy of S U{2) corresponding to the /-th edge. 

It is clear that the operator (2.4) depends on the triangulation T{e). It turns out that the nontrivial action of on a cylindrical 
function fy corresponds to v(A) ny 0, which motivates us to triangulate Z adapted to y [9]. We denote the triangulation adapted 
to y by T{y). The assignment T{y) which we choose is as follows; 

(1) If V e y(y) is non-planar with valence bigger than two, a tetrahedron A is assigned to v so that v(A) = v and the three edges 
starting from v(A) are chosen as the starting segments (s/, sj, sg) of a triple of distinct edges {ej, ej, eg) of y with indecent 
tangents at v. 

(2) If V is coplanar, we add one more edge e incident at v to y such that its tangent at v is transversal to the tangents of edges 
of y at V. The new graph is denoted by y'. The three edges of A with base-point v(A) = v come from the starting segment 
of the new edge e and the other two edges of y with independent tangents incident at v. 

Then the action of the regularized Euclidean Hamiltonian constraint operator (2.4) on fy reduces to 

H^yiN) ■ /, H^y.yfN) ■ fy = 2 N{v)^ I, ' frI, ' /r > (2-7) 

^ veV(r) ^ v(A)=v veV(r) 

where E(v) is the number of non-planar triples of edges of y or y' at v, N[, 

Hy Yj f^^^tJ'(/i„„(A)/z.^(A)V/z;j(A)) ■ (2.8) 

.■(A)=v 

The limit e —» 0 can be taken in the Rovelli-Smolin topology. The label T for the triangulation T{y) can be dropped off since the 
final limit operator is independent of e. 


2.2 The action of Hy{N) on a trivalent non-planar vertex of y 

The action of the Hamiltonian constraint (2.7) is local in the sense that it is a sum over independent vertices. Therefore, we can 
concentrate on its action on a single vertex. Given a spin network state T ^4^A) on a graph y, we consider a trivalent non-planar 
vertex v e y(y) at which three edges ei,e 2 ,63 incident. The terms in //(■4) directly associated to v can be represented by ' 


rffA)-.^(iy) 

yj,i 






(2.9) 


^For calculational convenience, the spin network states considered here are not normalized, although their intertwiners are normalized. They can be normal¬ 
ized by multiplying n/=i 7/ + 1. 
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M=0 


where denotes the intertwiner associated to v. The summation in the expression of in (2.8) 

is over only one tetrahedron A adapted to y at v. We will omit the notation A. Then the action of i/f on T'’ ^^(A) can be explicitly 

y-j,i 

written as 


■ r ,U) = ■ T' ,U) = ^ ■ T’ ,iA) 
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I TTC^ I T_T‘ 
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)-T'\XA), 

/ y,.hi 


( 2 . 10 ) 


where = [^i/2(/!)]'^b' Note that applying to T'' ^^(A) involves the actions of the holonomy and the volume operators. 


y-j.i 


Before calculating the action (2.10), let us recall the graphical method introduced in the first paper of the series [17]. The 
3y-symbol is represented by 


71 72 73 I ^ 

nil ffi2 ^3 


( 2 . 11 ) 


The orientation of the node is meant the cyclic order of the lines. A clockwise orientation is denoted by a sign and an anti¬ 
clockwise orientation by a “-h” sign. Rotation of the diagram does not change the cyclic order of lines, and the angles between 
two lines as well as their lengths at a node have no significance. Summation over a magnetic quantum number m is graphically 
represented by joining the free ends of the corresponding lines. The intertwiner associated to v in (2.9) is represented in graphical 
formula as (see Appendix A in [17]) 




“ \ Jljl.h )n 
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mi 7712 7773 


7771 '^'^2 


mi 


^ -.=7W = 0 - ,3 


where du = 2]^ + 1, and in the last two steps we have used the following two identities ((A.43) and (A.51) in [17]) 


( 2 . 12 ) 
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here a line with an arrow on it and a line without arrow denote respectively 


^U) _ — (_]\J 
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- (-!)■' 


6Z, = S, 


m,m' — rn- 


(2.14) 


The matrix element [7ij(he)]"'^ in spin j representation of the holonomy hg = hAA) is denoted by a blue line with an arrow on it 
as 


M/r,)]” - ra. 


(2.15) 


The orientation of the arrow is from its row index m to its column index n. Then the matrix element [njihg *)]"„, can be represented 
by 


[nj(h:^)]\ - [njihg-.)r„ = , [ 7 :j(hg)r\,q!; = n 


-m = n- 


(2.16) 


Notice that the holonomy acts as multiplication. Hence its action on the spin network states will involve the Clebsch-Gordan 
series, which can be represented graphically by 
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With these preparations, the part of y associated at v and the spin network state T'’ in Eq. (2.9) can be represented respec- 

7,J,i 

tively by 


ei 






(2.19) 


where we have omitted the the remaining notations ni,n 2 ,nj in the algebraic form (2.9). By introducing three pseudo-vertices 
vi,I - 1,2,3, we subdivide e/ into two parts s/ and Ij such that e/ = s/ o Ij and s/ = s/(A) matching the triangulation T{y). Then 

T'' „4A) in Eq. (2.9) becomes 

y-jJ 




( 2 . 20 ) 


where intertwiners associated to v/. Hence the original graph and the corresponding spin network state in 

(2.19) become 



We can also single out the corresponding part of T'’ -^JA) which only involves the holonomies /ij,, and denote it by T \XA) (the 

yj.‘ yj,' 

notation s denotes the segments s/), i.e.. 


T'XiA) = (U) 
y-j,i 


m\m2ins 





( 2 . 22 ) 


Now let us to calculate the action of the first term in (2.10) on T'X(A) via the graphical method. The action of [h^X a 

on T'XXA.) can be represented by 

yj,‘ 
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where = Ji ± 5 , and we have adjusted the coefficients sucfi tfiat tfie intertwiner /' at v is normalized as 


‘.-'f’y t 


= (-!)"'■ Jd 
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= Jd 


— ^Jda2da, 


(2.24) 


■ a, = J, a,=j, J = 0 


Recall tfiat the volume operator (2.5) vanishes coplanar vertices. Hence it has non-trivial action only at v, not v/. Its action in 
(2.23) reads 


y 






4^2 


.h . 


(2.25) 


where the operator qfj^h corresponds to the edges ii, S 2 , 53 with spins respectively. Notice that quK changes only the 

intermediate momenta «/,■■■ ,aK-i between ji and Jk of the intertwiner. In our case, the operator qj\j 2 i 2 . and hence V change 
02 = ;i,a3 = 73 into 


«2 = /i ± 
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iji-lju for/j=;i-i 
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03 = 73■ 


Hence we have 
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is 
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(2.27) 


where lo^, a'^} and |a 2 , 03) are the simplified forms of the intertwiners corresponding to two different chosen intermediate momenta 
for the fixed anglular momenta 5 ,/p 72,73 and tfie resulting 7 = 0. For given values of four spins 5 , 7 ^, 72 , 73 , tfiere are two 
allowed combinations of intermediate momenta (2.26). Hence the corresponding intertwiner space associated to v has dimension 
2. Furthermore, the volume operator is automatically diagonal on the 2-dimensional intertwiner space. This fact was pointed out 
in [21, 11], which will also be presented in Appendix A. Therefore we obtain 
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Hence we have 






(2.31) 


where in the fourth step we have changed the orientation of two arrows with spin by the rule (B.7), and then used the rule 
(B.8) to remove three arrows with the same orientation joint with a 3y-symbol, and we also used the rule ((A.54) in [17]) 

(2.32) 

to remove a loop in the fifth step. Thus the action of ^ on r'’’^^(A) is given by 
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(2.33) 


where in the second step we have used (2.16), and in the last step we have used the two identities (see Appendix B for proof) 

^ _(_l)t.+/,+)(_l)23+2yi pi 5 /ijpi 2 {'Ipi,/, (2.34) 
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Eq. (2.33) enables us to write directly down the results for (/, J, K) 6 {(2,3,1), (3,1,2), (1,2,3)) as 



fpfj 



(2.36) 


where 

H{f„ fj, Jk) = Yj 
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(2.37) 


Taking account of the identity (-1)^^^’^' = (-1)^*, the action of on r''’L(A) can be explicitly written down as 
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3 The inverse volume operator 

The Hamiltonian of a massless scalar field reads 


H^(N) = i J 
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(^) = ^ [HkinAN) + //der,0(A^)] , (3.1) 
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where n is the momentum conjugate to (p. In order to quantize H^(N) in framework of the scalar field coupled to gravity in LQG, 
the term can be regularized as [ 10 ] 
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where we have inserted 1 = [det(e[j)]^/ -^det(q')j in the second step, used e‘^(x) = ^{A'J(x),V{x,e)) and absorbed y( 2 c, e) : = 

^d&t{q){x) in the denominator into the Poisson bracket in the last step. Again we introduce a triangulation Tiy) of S adapted 
to a graph y. For a given tetrahedron A and its edge i/(A) sj, by the identity 


r dk{A[,(2c), y(v,e)2) = 2tr(T,/!/{/i/,y(v,e)h)+ o(e^), /i/=/ i^ka) , 
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Eq. (3.2) can be reduced to 
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Replacing n by -i%K6l6(f>, Poisson brackets by commutators times lf{iK), and substituting y —> y, H\^in,<piN) can be quantized as 
HMn,^(N)y^ ■ 632 '^ 4 ^ 4 g 6 Yj N{v'')X(v")X(v"') x,(v'',v''')xPv'',v)x,(v'',v') 
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jinjMnjjv=v' 


where X{v) j [2fR(v) + Xi{v)] is the sum of left and right invariant vector fields acting on the point holomomies U{v) defined 
in [ 22 ], and 

y,{v) tr(r,/!/[/r 7 \ y!]) = -tr(r,/!/yj / 17 ') (3.6) 

is e-independent for sufficiently small e. For sufficiently small e, the three characteristic functions in (3.5) vanish unless v = v' = 
y" = y'". Taking the limit e ^ 0 yields 



where we have used e, ,vte/mn = 3!d[,i5md^n. The operator y"’aitv is the quantum version of — -J =— = up to a constant, and 

J U ’ 6^ •Ydet(.?)(A) 

thus it is called the inverse volume operator. By introducing the manifestly gauge invariant operators [15] 

§/y(v):-d,/i^'(y)ij^>), (3.8) 
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the inverse volume operator V"'ait,v c™ be represented in terms of qij{v) as 


V-\\uv-fy^ ^ e‘^’^e^’^^qiL{v)qjM{v)qKN(v)- fy. 
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For the convenience of graphical calculus, one usually introduces the spherical tensors (or the irreducible tensor operators) 

(/r = 0, +1) corresponding to r, {i - 1,2,3) by 
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where we have used the following identity in the second step (see Appendix B.l in [17] for proof) 
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and the following identities in the third step 

(r,)\ = (v)^Cf;;', (3.13) 

here the overline denotes complex conjugation. 

In what follows, we consider the action of y"'aity on T'’’iXA) at a trivalent non-planar vertex v. Notice that the intertwiner 

_ ’ yJ,‘ _ 

space associated to v, which will be acted by y^'ait.v, is of one dimension. Hence the gauge-invariant operators §/y(v) and y^'ait.v 

take eigenvalues on the orthonormal spin network states of T' ’7^(A), 
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(3.17) 


In order to obtain the eigenvalues Qu, we need to calculate the action of ^ 2 ^(v) on T (A) or T rXA). In what follows, 

y,j,i yJ,‘ 

we only display the derivation of the two quantities Qu and Qu, and the remaining components of Qu can be written down 
similarly. 

Now let us consider the action of *^ 2 ^(y) on T''’A(A). Notice that the spherical tensors Tu can be represented by (see Appendix 

y,j,‘ 

Ain [17]) 


Hence we have 


i^j(^,)]^B = 2 V27'(27 + 1)(2; + 2) 


r + j 


'-/2j{2j+\){2j + 2) 



(3.18) 
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.V6 








(3.19) 


where in the second step we have used the result of Eq. (2.31), and in the fourth step used the identity (see Appendix B for proof) 




(3.20) 


Taking account of 


2 /;i(;i + 1) ^ jl, j\-j\ + \ 




Eq. (3.19) can be reduced to 


<2^(y) 



V;i(;i +1) I k _yi 


2ji + 1 


d,i- 



(3.21) 


where 




V'^B [V{j\ = A - 1/2,72,73)]^ . 


(3.22) 


The intertwiner in Eq. (3.21) is normalized because of 


dn (-1)"'^ 


h 
. ii' 


-M - 


yjdjj dj^ J 


Jl 

1 

32 

h 

K 





■ 02 = jl 

■ fls = js 

■ J = 0 


(3.23) 


Eq. (3.21) implies that eAv) changes neither the graph nor the spins of r'’l^(A) (2.22). But it does change the intertwiner 
into i[, = ;-2=0’‘'2=31.03=23 associated to v. Hence we obtain 

' 31 , 32,23 '' 31 , 1 , 32,33 


gii = - ((i^(v)r’'’l’r”, '5^(v)r'’'l'r”l - -iAMmi (vl _ yM _. V2i(A + / i _ i \ j,v,3,n, 

V ^ 7,3,1 1 y,j,i 2/1+1 \ 7,3,1' 2/| + 1 \ 7,3,1' 


. VTlUTTl) / I 1 


2Ji + 

30- f n d//H(/i.,)r’'’tf™(A)r'’’'^r”(A) 

( 2 ;i + 1)2 \ Jsj/(2)3 I 3 7,3,1' r.3,,' 


2Ji + 




;i(;i + 1) 

JiUi + 1 ) (,A 

'(2;i + l)2 


/yJ - y2 ') tr(/'^“°’ 


-0;a2-3l,03-33 _ y-0;02-31,03-33 I 
1 , 32,33 J\,ljl,h 
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(3.24) 


;i(;i + i)/ i _ /.y=0;a, 

yiA 

iy'^ - yM 

lA ^ \Bj ’ 


(2;i + 1)2 

j\{j\ + 1) i 


=yi,^'3=73 , ^J^^',a2=jv,a2-h\ 
•h 714 J 2 J 3 ) 


(2;i + 1)2 

where tr() denotes contracting magnetic quantum numbers, we have integrated holonomies to give the contraction of the inter¬ 
twiner with its complex conjugate in the fourth step, used the fact that the intertwiner is real in the fifth step, and the intertwiner 
is normalized in the last step. 

Similarly, the action of e'^(v) on ^'’’^^(A) yields 

yj,' 





, y/jiUi + 1) (i i 


2j2 + 1 


- yM 

yiA 



■i "‘yw (4 - 4 ) (-1)2-^2..A ^ + 1)(2;2 + 1) 


73 ]2 a 

ji h 



(3.25) 


where 

^2A =72 + 1 /2, ii , 73)] I , [y (72 = 72-1 /2, ji , 73)] ? , 

and in the second step we have used the following identity (see Appendix B for proof) 



1 + 72+73 


l,(2» + l)lf “ 


7 i 72 


a '' 


(3.26) 


(3.27) 


Hence ’^ 2 ^(v) changes i,, s i'.^ ^7 Jjjj-q composition of i" s associated to v. Finally, we have 


Q 


12 


\ ‘ rj.' 2 y,j,i 

__y[Mh^(y{ j \ V 72(72 + 1) / i ^y.+7: 

“ 271 + 1 \ W 272 + 1 \ 2+ ^ 


+ 72+73 


X y V(2a + 1)(272 + 1), 1 . . 

I i 7 i 72 

a ^ 

yJjiUi + 1)72(72 + 1) /, i 


73 72 a 


tr 


V 7i4 


=0; 6 ( 2 = 7i >< 23 = 73 ■/=0; a 2 ~a,as= 73 


72 J 3 


• I 


7l>l'72^3 


(271 + 1)(272 + 1) 

V 7 i( 7 i + 1)72(72 + 1) / i , A\(^A t/5 


(yL - 4)(4 - Z V(2«+ 1 x 272 +i)r; 

a ^ 

, 21 , . 1X22, . 1) (A - A)(A - 4)(-lV'*J-'>V(22, + l)(22,+ l){f 

V7i(7i + 1)72(72 + 1) i ,A\/,A 


73 72 a 

1 71 72/"“' 2 ' 


73 72 7 i 

7 i 72 


(4 - 4)(4 - V( 27 i + 1)(272 + 1) 


(271 + 1)(272 + 1) 

X (--lyi+22+73+1 2[7i(7j + 1) + 72(72 + 1) - 73(73 + 1)] 

V27i(27i + 1)(27i + 2 ) 272(272 + 1)(272 + 2) 
71(71 + 1) + 72(72 + 1) - 73(73 + 1) l^A 


. (y2 _ y5 Vy5 - yM 
y lA '^iBjy2A '^2Bj ■ 


2(271 + 1)(272 + 1) 

Similarly we can write down the remaining components of Qu and thus the eigenvalue of V"'ait,v in (3.16). 


(3.28) 
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4 Summary 


In this paper, the graphical method is employed to compute explicitly the actions of the Euclidean Hamiltonian constraint and 
inverse volume operators on the spin network states with trivalent vertices. The rules of transforming graphs in our method 
simplify greatly the calculation. Every graphical reduction in our graphical calculus corresponds to an algebraic reduction in an 
unique and unambiguous way. This ensures the rigour of our calculation. 

The action of the Euclidean Hamiltonian on the spin network states with trivalent vertex v was shown in (2.38). 

yjj 

The difference between our result (2.36)) and (11.14) in [7] is the factor -1/2. Based on the gauge invariant operators qijiy), 

the inverse volume operator y"'aitv defined in (3.9) takes eigenvalues on the orthonormal spin network state r'’l’"°™(A). The 
_ ’ _ yj'.i 

eigenvalues of presented in (3.16), consist of the eigenvalues Qu of qij{v). The eigenvalues of y"*ait,v were also derived 

in [15] by the algebraic method. The difference between the value of Q\\ in (3.24) and (4.17) in [15] is a minus sign, while the 

difference between the value of Qn in (3.28) and (4.17) in [15] is the factor -[2 + (-l)^^''‘’^^^^]/3. 

In principle, the graphical calculation method can be applied to the general cases where the spin network states are defined 

on arbitrary valent vertices and the holonomies appeared in the two operators are expressed in arbitrary representation of the 

gauge group. However, for those general cases, the volume operator lacks the explicit matrix elements formula. This prevents us 

from doing further calculation. Eor the same reason, the matrix elements of the Lorentzian part contained in the full gravitational 

Hamiltonian constraint operator have not been explicitly written down even on the trivalent vertices [23]. 
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Appendix A The diagonalization of volume operator in 2-d intertwiner space 

We denote 



j'l 

32 

33 


“2 = ^ 

fl 

■ “3 = h 

32 

■ J = 0 

33 

“2 =+ J 

B 

II ’ 

1 

II . 

0 


|ai) = |l/2,/i,72,73;a2 = h - 1/2,03 = 73,-/ = 0) = yjlai + 1 V2a3 + 

lai) = |l/2,/i,72,73;a2 = /i + 1/2,03 = 73,-/ = 0) = V2a2^V2a7^ 

The matrix of the operator iqfj 2 h I*' '■1*^ above two states reads 

(■'' ^<«'il'§7;7273l«2>\ / 0 {ai\iqfj2h\Q:2}\^IO -ib 

\^a2\iqj>j2j2\ay) {a2\iqf^j2h\^2)] \{a2\iqj>j2h\^,) 0 ] ' \ib 0 / ' 

The eigenvalues and corresponding (normalized) eigenvectors of iqj[j 2 j, are given by 

Hence we obtain 


■^\iqj\j2j3 \ lo'i) = ^ ■^\iqj\j2j3\ lo'i) = ^ |e,)<e/|ai) = Vi^ Icri), 

1=1 i=l 

2 _ 2 

■^\iqj\j2j3 \ l®2) = ^ ■^\iqj\j2j3\ l®2) = ^ Vw ki><ei|a2> = Vi^ \ai) - 


1=1 1=1 
Now we derive the value of \b\. Using the matrix elements of {a'\q22>A\3) in [17], we have 

{ai\qf^j2j2\a2) ={02 - 1 = /i " I\qfj2j3\a2 = /i + 1) 


(A.l) 

(A.2) 

(A.3) 

(A.4) 

(A.5) 

(A.6) 
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.i(_l)^/.+>3(_l)/. + ^(/.4)x0';,;2)5X02,73)^ ^2j[(2j\+2)(2h + l)|j 1 I 


(_1)/i-5+aU3 /1 + 2 U /1 

I ^ 2 -^2 I I 1 


2 J2 J3 
73 h 


(_l)7i+5+;3 J 73 J2 7i + 2 I J7i 
11 h~ 2 72 / \ 


73 73 

73 73 

72 73 

73 72 


= - 4 (- 1 ) ^ (-1 )X{j\, jl) ^ ^( 72 , 73 ) ^ 

X J2;;(2;;+2)(2;3 + l)fJ 1 t ^ 

' I i 7 i 2 '' 


7i 


0 73 731 j73 jl j'l + 2 

72 


1 73 73/11 7 ;-i 


/7i 2 72 73l_/_,y;+i +23 f7i + 2 72 73 

^ ^ I 1 73 72 I ^ ^ II 73 72 


73 Jl 


With the following values of 6y-symbols 


1+1 


! rfU = 

0 73 73 I _ (-1)'^'^"‘ 


2;;(2;; + 1)2(2/i+2) 


1 73 73 j 273 + I 

''' • 12(/j + jl + jl + |)(/j + jl - jl + j)(j\ - jl + jl + j)(-j[ + jl + jl + 2 ) 


2 I — ('_'I v1+72+73+5 
1 j\-\ 72 / ^ ^ 

7i 2 ^ 73/ _ /I \7 'i+72+73+5 
1 73 72 ^ ^ 


2^(72,/i) 

2[72(72 + 1) + 73(73 + 1) - (/i - \) U \ + 5)] 


2f(72,73)'^2 


7'i + 5 72 73/ ^ ,_,y;+22+23+3 2[72(72 + 1) + 73(73 + 1) - (/i + 2)(7'i + 2)] 

1 73 72/ 2f(72,73)'/2 

where X{j\,ji) = 2ji{2j\ + \){2ji + 2)2ji(2ji + l)(2ji + 2), the matrix element formula (A.7) can be simplihed as 


I, 


1 , 


<fl2 - 1 = j'l - = fi + 2 ^ " 


. 3 . . 1 . .1 . . .1 

(J'l + jl + 73 + o)(/l + 72 - 73 + o)(/l - 72 + 73 + o)(“/l + 72 + 73 + ^ ) 


2 ' 2 

which implies that the absolute value of b in the matrix elements (A.3) takes the form 


3 ... 1 ... 1 ... 1 

(j'l + 72 + 73 + + 72 - 73 + “ 72 + 73 + 2^(~fi + 72 + 73 + 2 ^ 


Therefore, we have 


VlaJ = ^|/§/,22aI la,) = ^ la,) 


4 V 2 

4 V 2 


4 V 2 

3 ... 1 ... 1 ... 1 

(/i + 72 + 73 + 2^(j'i + 72 - 73 + 2^(j'i - jl + jl + 2^(~fi + 72 +73 + 2 ) 


|a/) 


= V(ll2,j'i, jl, jl ; ai = /j + 1 /2, 0,3 = 73 ) Icr,), / = 1,2, 
which reveals that the volume operator is diagonal in the 2-dimensional intertwiner space. 


(A.7) 


(A.8) 


, (A.9) 


(A. 10) 


(A.ll) 


Appendix B Proof of some graphical identities 


In the graphical calculus, one usually uses the following identity ((A.60) in [17]) to simplify graphs. 



(B.l) 
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The graph on the left-hand side of (B. 1) can be transformed to 



where we have used the rules of transforming graphs ((A.50)-(A.53) in [17]) 


-m m- 


-m m- 


.j. , = i-i')2f 


= (- 1 ) 


2j 



Graphically, Eq. (2.34) can be proved by 




j'' ’ 







(B.2) 


(B.3) 


(B.4) 


(B.5) 

(B.6) 

(B.7) 

(B.8) 


(B.9) 


(B.IO) 
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(B.ll) 


=-(-iy‘+^'>+n_iy3+;3+i O2 .^3 {1X1^2 h j}\ J\., 
Ul 2 ^ 2 ] I 2 Ji) 


^_(_iyi+/.+^(_iy3+2yU(i 2 yly) 2 V , . 

[J 2 J2 J2j Ua J 2 JsJ 

where we have used the rules (B.5)-(B.8), and (B.IO) in the third and fifth steps respectively. Similarly, Eq. (2.34) can be shown 
by 




— ^ ^J2 1 >^'3 y'l 1 l f 

l/i j's 2 J u'l 2 Al 

= (_iyi-/i+5(_iy2+22+5 i-fi 2 /ilyi 2 /il .yil ., 

{J 2 A A) u '3 }2 h] ■’>■>'^'1-’ 


(B.12) 


where we have used (B.IO) in the second and fourth steps, and used the fact that the allowed triple j'^, j\) satisfy the triangular 
condition in the fifth step. 

The identity (3.20) can be proved by 


-Li i 


Jlli J. 

yi 2 


- ji 

= (_l) 2 /,(_l) 2 xi y 


-M = (-l)22i + i 



= (-l) 2 /i + l I 2 I 2 1 = (-l) 2 /i + l M 2 A\ 

[ji Ji ji) ii| (1 ji 2) yy 

^(_ 1)2/,+ 1(_ 1)221 + 1 5 i| i|j_^ = (_l)2/+22,+2|;; 5 ^/|-|^ 


A i y p 
1 A Uj 




(B.13) 


where we have used (B.l) in fourth step, and used the fact that (-l)2/i+2/+i = l in the last step, since the allowed triple 5 ) 

satisfy the triangular condition. 

Eq. (3.27) can be obtained from 



u 


'3 = +; 


4^ = y(2fl+ !)(-!)' 

t J 


1 +/+72+ 0)72 73 O 

71 1 72 


1 +_/3 + 72 +« 1 ^ 


= J](2fl+l)(-l)i 

a 

= (-1)71+22+/ Y^(2a + 1) 

a 

where in the second step we have used the identity ((A.62) in [17]) 


7i 1 72 


“ l,(_l) 7 i+«+l(_l) 2 « 


{73 

72 

a \ 

i 


11 

7i 

h] 

a' 



mi , , 7714 

= y(2;6 + l)(-l)72-"73+/+;6|77 72 76 

/j3 jX ^ 1,74 73 75 

m 3 7112 


mi m4 

Vii uy 


(B.14) 


(B.15) 
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